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We have investigated several properties of rapidly rotating dynamic black holes generated by 
gravitational collapse of rotating relativistic stars. At present, numerical simulations of the binary 
black hole merger are able to produce a Kerr black hole of Jfinai /Af| Ilal up to = 0.91, of gravitational 
collapse from uniformly rotating stars up to Jfinai /A/fLai ~ 0.75, where Jfi na i is the total angular 
momentum and M nna i the total gravitational mass of the hole. We have succeeded in producing 
a dynamic black hole of spin Jfiiiai/A^Jnai ~ 0.95 through the collapse of differentially rotating 
relativistic stars. We have investigated those dynamic properties through diagnosing multipole 
moment of the horizon, and found the following two features. Firstly, two different definitions of the 
angular momentum of the hole, the approximated Killing vector approach and dipole moment of the 
current multipole approach, make no significant difference to our computational results. Secondly, 
dynamic hole approaches a Kerr by gravitational radiation within the order of a rotational period 
of an equilibrium star, although the dynamic hole at the very forming stage deviates quite far from 
a Kerr. We have also discussed a new phase of quasi-periodic waves in the gravitational waveform 
after the ringdown in terms of multipole moment of the dynamic hole. 

PACS numbers: 04.25. dg, 04.25. D-, 04.30.-w, 04.40.Dg, 97.10.Kc 



I. INTRODUCTION 

There are various mass ranges of black holes (BHs) in 
nature. Supermassive BHs exist in the center of most 
galaxies, and the typical mass range of this category is 
around 10 5 M Q - 1O 1O M . In spite of clear evidence of 
its existence, the actual formation scenario for the su- 
permassive BHs is still not certain [1]. There are sev- 
eral candidates for the intermediate mass range (around 
10 2 M Q - 1O 3 M ) of the BHs in the globular clusters @. 
At present the object of the intermediate mass range has 
not been directly found yet. Moreover, the standard for- 
mation scenario for such an object has to pass through 
the formation of stellar mass objects. The merger of stel- 
lar holes or compact objects, or collision and collapse of 
massive stars is a typical scenario for forming such inter- 
mediate mass range objects. There are also some other 
candidates for stellar mass range (around 3M Q - 50M Q ) 
of BHs in our galaxy. Binary coalescence of the stars or 
the collapse of the star of stellar mass range is a typical 
scenario for forming such objects. 

Nowadays, we can produce a dynamic BH by com- 
puter. We have two representative scenarios for form- 
ing a dynamic BH promptly. Here we neglect accretion, 
since the standard timescale for this process is consider- 
ably longer than the dynamical one. The first one is the 
merger of equal mass, binary BHs. Based on the cur- 
rent numerical simulations of the binary BHs, there may 
exist an upper spin limit for the newly formed BH. The 
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binary composed of non-spinning individual BHs leads 
to the final maximum spin of the newly formed BH of 
•/final /Mimd = 0-69 (Jfinai is the angular momentum and 
•Mfinai the gravitational mass of the newly formed BH) 0] , 
while the spinning individual BHs in arbitrary direction 
leads to a final BH spin of J nnal /M| nal = 0.91 0. More- 
over, test particle approximation in BH perturbation ap- 
proach including the superradiance effect leads to a final 
BH spin of Jfinai/M| na i = 0.9979 5]. In theory, we have 
the following discussion to support the existence of the 
upper spin limit of the newly formed BH. If the plunge 
phase of the binary BHs is characterised by the physi- 
cal quantities at a certain separation radius, namely the 
ISCO (innermost stable circular orbit) of a newly formed 
BH, then there may exist an upper limit to the spin of 
the newly formed BH because most likely there exist a 
radially unstable condition at the ISCO under which the 
binary begins to collide and form a new BH. 

The next one is gravitational collapse of a uniformly 
rotating relativistic star. In this scenario, the maximum 
spin of the BH exists by the following discussion. First, 
a star contracts itself to the mass shedding limit, con- 
serving the angular momentum of the system. Then, 
so far as the system contains sufficient angular momen- 
tum, the star evolves along the mass shedding sequence 
quasi-stationary, releasing the mass and angular momen- 
tum. Once the star reaches the critical onset of collapse 
because of relativistic gravitation, it begins to collapse 
6] . From the collapse of a uniformly rotating supermas- 
sive star, the final spin of a newly formed BH is around 
Jfinai/M| nal a 0.75 0. 

One of the primary observational missions for detect- 
ing gravitational waves in ground-based and space-based 
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TABLE I. Two different radially unstable rotating equilibrium supermassive stars for a BH formation 



Model 


Rp/ R£ 


niaxb 

Po — 


M c 




J/M 2 2. 


M/Rl 


I 


0.450 


1.56 x 10~ 5 


4.88 


0.108 


0.99 


2.56 x 10" 2 


II 


0.425 


1.56 x 10" 5 


5.07 


0.118 


1.03 


2.63 x 10 -2 



a Ratio of the polar proper radius to the equatorial proper radius 
b Maximum rest mass density 
c Gravitational mass 

d Ratio of the rotational kinetic energy to the gravitational binding energy 
° J: Total angular momentum 
f R: Circumferential radius 
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FIG. 1. The gravitational mass and angular momentum 
of the dynamic BH through gravitational collapse for model 
I. The angular momentum is computed by the approximate 
Killing vector (solid line), while by dipole moment of the 
imaginary part of the Weyl scalar vE^ (dashed line) on the 
apparent horizon. The gravitational mass is computed by 
using the first law of BH thermodynamics for both cases. 
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FIG. 2. Same as Fig. Q] but for model II. 



interferometers is to investigate a various mass range of 
BHs and compact objects [Sj. Combining the global net- 
work of gravitational wave detectors, we are in these 
decades able to extract fruitful features of BHs in the 
frequency band of 10 _4 Hz - 10 3 Hz. Potential sources of 
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FIG. 3. Euclidean norm of the normalised Hamiltonian con- 
straint and the normalised Momentum constraints throughout 
the evolution for models I (solid line) and II (dashed line). 



high signal to noise events in this frequency range are 
quasi-periodic waves arising from nonaxisymmetric bars 
in collapsing relativistic stars and from the inspiral of bi- 
nary BHs for example (e.g. Ref. §)■ In addition, a non- 
spherical collapse of a rotating relativistic star to a BH 
potentially generates a significant amount of burst waves 
and quasi- normal ringing waves (e.g. Ref. [13]). In this 
paper we trace the collapse of relativistic stars through 
numerical simulations to investigate some of these possi- 
bilities. 

Here we relax the condition of uniformly rotating pro- 
file in the equilibrium to produce a highly spinning dy- 
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namic BH. Differential rotation profile of the star enables 
us to impose large amount of angular momentum in the 
system, since it relaxes the restriction to the angular ve- 
locity at the equatorial radius, which comes from the 
limitation of the mass-shedding. According to the above 
idea, Saijo and Hawke [Hi have succeeded in producing 
a dynamic BH of spin Jfi na i/Af| nal s» 0.98. Here we focus 
on the BH configuration in this paper by using multipole 
moment of the curvatures on the apparent horizon. We 
try to answer the following questions. Can we extract 
precisely the mass and angular momentum of a dynami- 
cal BH by using multipole moment of the curvatures on 
the horizon? Can a newly formed BH be represented 
as a stationary Kerr BH at several dynamical times af- 
ter the BH formation? Is it useful to use multipole mo- 
ment of a dynamic BH to extract some properties of a 
BH, and to find a cause of quasi-periodic gravitational 
waves after the ringdown, for example? To answer these 
questions, three spatial dimensional general relativistic 
hydrodynamics is necessary. 

The content of this paper is as follows. In Sec. [Til 
we briefly explain the general relativistic hydrodynam- 
ics, especially the numerical tools we use to understand 
the property of a dynamic BH. In Sec. IIII1 we introduce 
our findings of a dynamic BH, focusing on its configu- 
ration. Section IIVI is devoted to the summary of this 
paper. Throughout this paper, we use the geometrized 
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units with G = c = 1 and adopt Cartesian coordinates 
(x,y, z) with the coordinate time t. Note that Greek in- 
dex takes (t, x, y, z), while Latin one takes (x, y, z). 



II. BASIC TOOLS IN NUMERICAL 
RELATIVITY 

In this section, we briefly describe three-dimensional 
relativistic hydrodynamics in full general relativity. We 
also explain our techniques for investigating outgoing 
gravitational waves from the sources and a dynamic 
horizon configuration (see Ref. 11| and references cited 
therein) . 



A. The gravitational field equations 

We define a spatial projection tensor = g^ v +n^n v : 
where g^ v is the spacetime metric, = (l/a,~f3 l /a) 
the unit normal to the spatial hypersurface, and where a 
and /3 l are the lapse and shift. 

We evolve the spacetime with the 17 spacetime associ- 
ated variables (</>, K, 7^, Ay, T l ), where is the con- 
formal factor, Kij the extrinsic curvature, 7^ the con- 
formally related spatial 3-metric, Ay the conformally re- 
lated trace-free extrinsic curvature, and T l the conformal 
connection function. The evolution equations are 
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(2.1) 
(2.2) 
(2.3) 
(2.4) 

(2.5) 



where Cp denotes the Lie derivative along the shift /3 1 , 
Pn = T^n^n", Si = T^mPhX and TF the trace-free 
part of the tensor. This set of equations for solving the 
Einstein's field equations numerically is usually called the 
BSSN formalism. As for gauge conditions, we choose the 
generalised hyperbolic if-driver pjffor the lapse, and 
the generalised hyperbolic T-driver [13| for the shift. 



B. The matter equations 

We assume a perfect fluid for describing a relativistic 
star as 

= p(l + £+^u»u" + Pg^, (2.6) 

where p is the rest-mass density, e the specific internal 
energy, P the pressure, and it^ the four-velocity. 

Energy momentum conservation V M T MI/ = together 
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with a continuity equation, leads to the flux conservative and the relativistic Euler equations as \\A 
form of the relativistic continuity, the relativistic energy 



1 d 

/=gdt 
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i—g dx 



-{J—gF) = S\ (2.7) 



where the state vector U, the flux vectors F 1 , and the 
source vectors <S l are 



U= [D,S,,T} rj 



Pv l 



0,T' 



(2.8) 
(2.9) 

(2.10) 



and r* - fc is a Christoffel symbols. Note that (p,Vi,e) are 
the physical variables of the above equations, and the 
flux conserved quantities D, Si, r are 



D = P W 7 (2.11) 

W = au t , (2.12) 

S t = P hW 2 v l: (2.13) 

E = phW 2 - P., (2.14) 

t = E — D, (2.15) 



where Vi — Ui/W, h = 1 + e + P/p is the specific en- 
thalpy. In the Newtonian limit, the above three physical 
variables coincide with the rest mass density, the flux 
density of the rest mass, and the energy of a unit volume 
of the fluid. In order to solve the set of equations, we 
have to impose an additional condition among the ther- 
modynamical quantities, namely the equation of state. 
We adopt a T-law equation of state in the form 

P=(r-l)pe, (2.16) 

where V is the adiabatic index which we set to 4/3 in this 
paper, representing a supermassive star (the pressure is 
dominated by radiation). 

C. Gravitational waveforms 

We monitor the Weyl scalar ^4 in Newman-Penrose 
formalism for investigating outgoing gravitational waves 
as 

*4 = C^ Xa k»m v k x m\ (2.17) 

where C^uXa is the Weyl tensor, A: M the ingoing null vec- 
tor, m M and mf 1 are the orthogonal spatial-null vectors of 
the four complex null tetrad (/ M , fc M , m M , m M ). The Weyl 
scalar ^4 represents the outgoing gravitational waves at 
infinity 

94, = h+-ih x , (2.18) 



where h + and h x are the two polarisation modes 
(transverse-traceless condition) of the perturbed met- 
ric from flat spacetime in spherical coordinate, and q is 
the time derivative of the quantity q. The Weyl scalar 
^4 roughly represents the outgoing gravitational waves, 
ignoring the radiation scattered back by the curvature 
when locating the observer far from the source. Therefore 
we trace the Weyl scalar ^4 to understand key features 
of gravitational waves emitted from this system. 

D. Horizon configuration 

Here we introduce a useful idea to diagnose the hori- 
zon locally in dynamical spacetime. It is the dynamical 
horizon defined as the outermost trapped tube which is 
composed of the apparent horizon in our case [15j . First 
we have to define the angular momentum from the hori- 
zon configuration. One way to determine an angular mo- 
mentum of the dynamic BH is (see e.g. section III.B of 
Ref. [H) 

Jbh = ~ I Kr,R»<p v ds, (2.19) 
SlT Js R 

where is the outward directed spacelike normal to 
the horizon in the spacelike slice, and <p a is a rotational 
vector field on the horizon. This quantity is interpreted 
as the Komar angular momentum when (p a is a rotational 
Killing vector on the horizon. The code and method used 
to compute numerically an approximate Killing vector, 
should it exist, is described in Ref. (T7j . 

Another way of defining an angular momentum, which 
is dipole moment of the scalar curvature is 

where ^2 is Weyl scalar, Y l m the spherical harmonics, A 
the area of the horizon and u the polar angle of the hori- 
zon configuration. Note that the computations from two 
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FIG. 4. Gravitational waveforms for model I. We moni- 
tor the real part of the Weyl scalar ^4, which represents 
the plus mode of outgoing waves at null infinity. The ob- 
server is located in the rr-direction of the equatorial plane 
at x = 65.52M, 98.29M, 131. 05M from the top panel, re- 
spectively. Note that the apparent horizon appears in the 
hypersurface after t — 670. 97M. 
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FIG. 5. Same as Fig.0]but adjustment of time shift, plotted 
in the same panel. Solid, dashed and dash-dotted line denotes 
the waveform detected in the ^-direction of the equatorial 
plane at x = 65.52M, 98.29M, 131. 05M. We use the following 
adjustment of time shift t a dj = t + ^farthest observer — Xobs, 
assuming that gravitational waves propagate with the speed 
of light in flat spacetime. 
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FIG. 6. Same as Fig. [4] but for model II. The observer 
is located in the ^-direction of the equatorial plane at x = 
63.15A1", 94.72M, 126.30M from the top panel, respectively. 
Note that the apparent horizon appears after t — 770. 40M. 
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FIG. 7. Same as Fig.|6]but adjustment of time shift, plotted 
in the same panel. Solid, dashed and dash-dotted line denotes 
the waveform detected in the ^-direction of the equatorial 
plane at x = 63.15M, 94.72A4", 126.30M. We use the following 

adjustment Of time Shift t adj = t + ^farthest observer - Zobs, 

assuming that gravitational waves propagate with the speed 
of light in flat spacetime. 
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FIG. 8. The m = 1 diagnostics of the rest mass density along 
the equatorial ring for model I. We measure the diagnostics 
in the equatorial plane at the radius r = 1.024M, 2.048M, 
4.095A/, 10.24M from the top panel, respectively. 



different definitions of the angular momentum coincide 
with each other in an axisymmetric apparent horizon. 

The dynamical horizon mass of the hole Mbh is com- 
puted once we have extracted the angular momentum of 
the hole by the following relation 



Mi 



1 



BH 



2-Rbh 



-Rbh + 4J BH . 



(2.20) 



Note that i?BH is the area radius of the hole. 

Axisymmetric isolated horizons are represented by two 
types of multipole moment of the scalar curvatures on the 
apparent horizon as [l8[ 



L t = 
h = 



1. 



^ 2 (u)Y l (u)dA, (2.21) 
2 K{u)Y l (u)dA, (2.22) 



where is the imaginary part of the Weyl scalar ^2 

*2 = \c^ p (l»n v l x nP - l^i v m x mP), 

which represents the gravitational monopole at large ra- 
dius, and 2 1Z the Ricci scalar. The quantities L\ and 
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FIG. 9. Same as Fig. [S] but for the m = 2 diagnostics. 



correspond to the mass and current /-pole moment de- 
fined in axisymmetric hole as 



Jl 
Mi 



4?r 

21 + 1 4tt 1 



4tt M BH i?BH 



21 + 1 



2tt 



Ji. 



(2.23) 
(2.24) 



This method has two disadvantage when computing mul- 
tipole moment numerically. One is that the quantities are 
gauge dependent, and the other is that they become less 
accurate in the fixed grid of finite differencing when com- 
puting higher /-pole moment. In order to avoid the above 
two issues, we use a different method for computing Z-pole 
moment of the curvatures by introducing averaged quan- 
tities on the trapped surface. We introduce the following 
'V-pole moment of 2 1Z and as [l|| 

Hn( 2 n) = ((( 2 K)~ 2 KT), (2.25) 

M$** 2 ) = («3*2) - 9f*2) n >, (2.26) 

where the bracket of a physical quantity (Q) represents 
the averaged quantity on the BH horizon 

{Q) = \J Q dA. (2.27) 

The definition of n-pole moment is a general extension of 
defining the variance of quantities. The relations between 
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FIG. 10. Same as Fig. [5] but for model II. We measure the 
diagnostics in the equatorial plane at the radius r — 0.987M, 
1.973M, 3.947M, 9.867M from the top panel, respectively. 



the quantities /z„( 2 7?.) and /^(G^), and the mass and 
current ^-pole in axisymmetric spacetime are 



^Cn) = ^l-2j2liY l o(u)) ), (2.28) 
Mn(^*2) = ( (jrLtY l (u) ) }. (2.2')! 




\2=0 



Therefore /j, n ( 2 lZ) corresponds to the summation of all 
mass ^-poles, and /^(Ss^) ai l current Z-poles. Since we 
impose planar symmetry across the equatorial plane, the 
current odd ^-poles vanish entirely. The disadvantage 
of using the quantities //„( 2 7£) and /^(O^), however, 
is that it is quite difficult to understand their physical 
situation. Therefore, we introduce the nondimensional 
quantities for 2 1Z and 3\I/2 as 



l n = 



8tt, 



A 

from a computational viewpoint. 



(2.30) 
(2.31) 
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FIG. 11. Same as Fig. UUI but for the m = 2 diagnostics. 



III. NUMERICAL RESULTS 

Here we focus on the properties of a dynamic BH 
through gravitational collapse of a supermassive star. We 
choose the equilibrium star radially unstable for evolu- 
tion in order to focus on the BH formation fill ] . 

First we compute the gravitational mass and angular 
momentum of the hole with two different definitions for 
the angular momentum in Figs. [T] and [2j One definition 
for computing the angular momentum is to use the ap- 
proximated Killing vector, while the other is to use dipole 
moment of the imaginary part of Weyl scalar ^2- When 
the dynamical system is axisymmetric, the computations 
of the angular momentum by two different definitions co- 
incides with each other. We find a clear agreement of the 
gravitational mass and angular momentum between two 
different definitions in Fig. [1] for model I and in Fig. [5] for 
model II. The results also tell us that our gravitational 
collapse is nearly the same as an axisymmetric dynamics 
for both models. 

Next we monitor the Hamiltonian and Momentum con- 
strains through gravitational collapse for monitoring the 
accuracy of our dynamics in Fig. [3] These checks are nec- 
essary because we do not solve these constraints through 
time integration of the Einstein's field equations. We di- 
agnose the same quantities as before |11| . the Euclidean 
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FIG. 12. Multipole moment of the Ricci scalar 7? of the dy- 
namic BH for model I through evolution. Solid line represent 
the one from our dynamic BH, while dashed line is the one 
from a Kerr using the gravitational mass and angular mo- 
mentum of the dynamic BH. We also enlarge the region of 
quasi-stationary stage at the right part in each panel. 
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FIG. 13. Same as Fig.Q/JJbut for the imaginary part of Weyl 
scalar $2 of the dynamic BH. Note that the odd Z-pole of 
are exactly zero when we impose planner symmetry across the 
equator. 



norm of HC and MC x>VjZ of all grid points outside the 
apparent horizon, through the normalisation of the max- 
imum rest mass density outside the horizon in the same 
hypersurface. The maximum violations from the con- 
straints of MC X and MC y are less than 1 x 10~ 2 for 
model I and 5 x 10~ 3 for model II. These facts tell us 
that our computational results are very accurate. They 
are roughly less than around 1% relative error at their 
worst. However, all the relative violation errors from the 
constraints for model I increase at the very late time of 
evolution (t m 900M ~ 1000M). We therefore stop our 
time integration at the time around t « lOOOAf to guar- 
antee roughly 1% relative error or less. 

We show gravitational waveforms using the Weyl scalar 
VfU through gravitational collapse in Figs. 0] and PoJ The 
Weyl scalar contains both outgoing waves and back- 
scattered waves by the curvature when we measure the 
quantity in finite radius from the centre. In order to 
focus on the outgoing waves, we monitor the waveform 
at three different locations, and investigate all of them. 
Since all three locations are considered as radiation zone 
of gravitational waves from the source, the gravitational 



field of all three locations is weak, and back-scattered 
waves only play a secondary role. The outgoing waves 
propagate towards spatial infinity as time goes on, the 
features of the outgoing waves can be seen in all three 
locations with positive time shift. We adjust the time 
axis of all three waveforms by assuming that gravitational 
waves propagate with the speed of light, and plot them 
in the same panel in Figs. [S] and [7] In fact, we use 
the following relation rj adj =t + a)f art host observer - ^obs to 
adjust the time. Although there is some difference in the 
magnitude of the amplitude, the global features look the 
same. Investigating all three waveforms, we find that the 
outgoing waves contain three features. The first is that 
there appears a burst wave as the collapse goes on. The 
second is that once the BH forms, there is a damping 
wave which corresponds to a characteristic oscillation of 
the dynamic BH. The third is that there is a continuous 
wave after the damping one. 

In order to identify the cause of continuous waves after 
the ringdown, we first investigate the azimuthal modes 
of the rest mass density. We introduce the following di- 
agnostics at certain radii of a ring in the equatorial plane 
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FIG. 14. Same as Fig. [121 but for model II. 
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FIG. 15. Same as Fig. [H but for model II. 



as 

Cm = 7T-FT- / P elmiP d ^ 

with a normalisation of D r ing(= Co), a mean density of 
the ring at certain radii in the equatorial plane. We in- 
vestigate m — 1 and m — 2 diagnostics at 4 different radii 
for models I and II in Figs. |8]-[TTJ Although the satura- 
tion amplitude for different radii is different for each m 
diagnostic, we find the following features. The azimuthal 
diagnostics begin to amplify efficiently after the apparent 
horizon has appeared in the hypersurface. This feature 
raises a question as to whether the amplification of the 
azimuthal diagnostics is directly connected with the con- 
figuration of the BH. The saturation amplitude of each 
m diagnostics is quite similar at the same radius. The 
saturation amplitude decreases as the radius becomes far 
from the BH. This feature suggests that the matter which 
is very close to the BH may play a key role for generating 
the quasi-periodic waveform after the ringdown. 

Next, we investigate the BH configuration to identify 
a possible cause of the continuous waves. We compute 
n-pole moment of the Ricci scalar 1Z and the imaginary 
part of the Weyl curvature $2 on the apparent horizon 
throughout the evolution. We also compute n-pole mo- 
ment of the same scalar curvatures using the configura- 



tion of a Kerr BH. We use the area of the horizon and the 
nondimensional Kerr parameter Jbh/^bh f° r computing 
n-pole moment. Then, we compare each n-pole in both 
dimensional definition with the BH area in Figs. [T21-IT51 
and nondimensional one in Figs. \TE\ - Q33 We find the 
following two features. The first is that after t lOOAf 
time from the formation of the dynamic BH Q, it is de- 
scribed as a Kerr BH is within a relative error of several 
percents. If we take the dimension of n-poles into ac- 
count through the area of the BH, n-pole of the dynamic 
BH approaches the one of a Kerr in Figs. [T^] and [T3] af- 
ter t = 750M, and in Figs. [H and [TS] after t = 850M. 
Therefore the BH configuration becomes nearly the same 
as a Kerr after t w 100M from the BH formation. This 
statement suggests that the cause of continuous waves 
may be related to the matter instability, since the BH 
configuration is nearly the same as a Kerr. The other 
is that the odd n-pole moment has large deviation from 
that of a Kerr. This may be accepted fact as the nonax- 
isymmetric configuration of the rotating BH, as it traces 
the violent phenomenon at the BH formation. One cau- 
tion from this feature is that the BH mass and angular 



1 We define the formation time of the BH as the first existent one 
of the apparent horizon in our hypersurface. 



10 



1 I I rq 




t/M 

FIG. 16. Multipole moment of the Ricci scalar R of the dy- 
namic BH for model I through evolution. Solid line represent 
the one from our dynamic BH, while dashed line is the one 
from a Kerr using the gravitational mass and angular momen- 
tum of the dynamic BH. 




t/M 



FIG. 17. Same as Fig.[l6]but for the imaginary part of Weyl 
scalar \l/2 of the dynamic BH. Note that the odd Z-pole of 
are exactly zero when we impose planner symmetry across the 
equator. 



momentum are settled down after t ps 100M from the 
BH formation time (See Figs. Q]and[2]), * na * the is 
almost regarded as a Kerr. For example, the half-life 
period of the BH oscillation r(= l/3u; qnm ) is 13. 6M, 
since the quasinormal mode frequency of a Kerr BH of 
a/M = 0.98 is Mu qnm = 0.422 + 0.0735i for I = 2, m = 
[20j . The fact leads to the conclusion that we cannot 
extract the "stationary" mass and angular momentum 
of the dynamic BH by quasinormal ringing in principle. 
Those ringing waves represent vibration of a transient 
dynamic BH, not a "stationary" one. 

IV. CONCLUSIONS 

We investigate the formation of the dynamic BH 
through gravitational collapse by means of three- 
dimensional hydrodynamic simulations in general rela- 
tivity. We particularly focus on the configuration of a 
dynamic BH and find the following two features. 

We investigate two different definitions for the angular 
momentum of a dynamic BH in order to check the va- 



lidity of the approximated Killing vector approach. We 
compare two results from two different definitions for the 
angular momentum and find that we are able to extract 
precisely the BH mass and angular momentum even if 
we use the approximated Killing vector. The fact also 
indicates that our cases of gravitational collapse are very 
close to axisymmetric. 

We also demonstrate the method to extract n-pole mo- 
ment of the dynamic BH precisely without using the ap- 
proximated Killing vector. This finding opens a new field 
of investigating the BH itself by extracting the properties 
of n-poles of the curvatures on the horizon. We compare 
the configuration of the dynamic BH with that of the 
Kerr, using multipole moment of the curvatures on the 
horizon. We find, as a result, that the quasistationary 
stage of the newly formed BH is approximately described 
by a Kerr BH. This does not mean, however, that the 
whole spacetime is approximately represented by a Kerr, 
since we only investigate the trapped surface of the hori- 
zon, just a local structure of the whole spacetime. 
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t/M 

FIG. 18. Same as Fig. [TBI but for model II. 
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t/M 

FIG. 19. Same as Fig. HH but for model II. 
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Appendix A: Multipole moment of the curvatures 
on the horizon in Kerr spacetime 



The 2-surface on the horizon of the Kerr metric in 
Boyer-Lindquist coordinate is given as 



■ sin 6 dip 



^[f-^ujdu 2 + f(u)dip% 



(Al) 





= r\ + 


a 2 cos 2 0, 




= M + 


VM 2 - a 2 


n 


= Kh 


-« 2 ) 1/2 , 




= a(rl 




u 


= cos 9 





/(«) 



1- P{l~u 2 )' 



(A2) 
(A3) 
(A4) 

(A5) 
(A6) 

(A7) 



The quantities 1Z and of the Kerr BH on the hori- 
zon are 



l n = - 



Stt d 2 |"(l + c 2 )(l-u 2 ) 
A du 2 



2tt d 2 
A du 2 



2(l + c 2 u 2 ) 
u{l + c 2 ) 2 
2c(l + c 2 u 2 ) 



where 



i-VT 



l + yr^a 2 ' 



(A8) 
(A9) 

(AlO) 



and a the nondimensional Kerr parameter. In order to 
compute multipole moment of the horizon only from the 
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nondimcnsional Kerr parameter and the horizon configu- vatures 2 1Z and • 

ration, we use the nondimcnsional quantities of the cur- Then, we can compute multipole moment of the Kerr 

horizon analytically as 



,1^ -15 - 70c 2 + 128c 4 + 70c 6 + 15c 8 3 ,„ ^ x4 arctanc 

^ = WT^) + Is (1 + c ) -i-' (A11) 



1125 + 5745c 2 + 10965c 4 - 1399c 6 + 6999c 8 + 6603c 10 + 2415c 12 + 315c 14 

2560(1 + c 2 ) 2 

+ 3^(1 + a 2 ) 4 (-25 + 14c 2 + y^^tanc 
512 c 

ln{ 2 Tl) = 1 ^(-294525 - 1755600c 2 - 4246935c 4 - 5382960c 6 + 4201406c 8 + 5703728c 10 + 6818482c 12 

+5017584c 14 + 2127279c 16 + 480480c 18 + 45045c 20 ) 

' 693 -(l + c 2 ) 4 (85 - 60c 2 + 22c 4 + 52c 6 + 13c 8 )^^, (A13) 



78848' 7 ' 7 c 

H^U) = - 75694Q8 | 1 + g 2 )4 ( 8437275 + 58760625c 2 + 175151130c 4 + 296175990c 6 + 325931705c 8 - 152055293c 10 

+104321660c 12 + 411065348c 14 + 548524189c 16 + 426449639c 18 + 207094602c 20 + 6 1 93 9 8 78c 22 
+10465455c 24 + 76 5 7 65c 26 ) 

:U(>:> (1 + c 2 ) 4 (-2435 + 1710c 2 - 1485c 4 - 572c 6 + 2067c 8 + 1326c 10 + 221c 12 ) ajctajlg ) ( A 14) 



7569408' 7 v ' c 

fi 6 ( 2 TZ) = X — — (-13841202375 - 111349888650c 2 - 396709663350c 4 - 835111292730c 6 

PbV ' 8921808896(1 + c 2 f y 

-1182872163330c 8 - 1230909094610c 10 + 906857870914c 12 + 1101124533086c 14 + 1994477265504c 16 

+3065485548066c 18 + 3312605946814c 20 + 2475795093330c 22 + 1279104290178c 24 + 449618015418c 26 

+102838870134c 28 + 13822058250c 30 + 829323495c 32 ) 

765765 

+ nr ^ 10 ^ 00 fl + c 2 ) 4 (l8075 - 11240c 2 + 15380c 4 - 920c 6 - 10302c 8 + 18216c 10 + 23252c 12 + 8664c 14 
8921808896 v ; v 

i c. arctanc 

+1083c 16 , A15) 

c 

, ~ . -15 + 170c 2 + 112c 4 + 70c 6 + 15c 8 3 , 2 . 4 arctanc 

= ^ra 2 ) + is (1 + ^ — a"' (A16) 

M4(3?* 2 ) = ,) — ^(-3465 - 36960c 2 + 619773c 4 + 663168c 6 + 1273910c 8 + 1306240c 10 + 985930c 12 + 515328c 14 

+178563c 16 + 36960c 18 + 3465c 20 ) 

693., „,, s arctanc < , . 

W 1+c) ^^- (A17) 

MeO^) = T^TT^^rr, ^(-271846575 - 4530776250c 2 - 35485039590c 4 + 755725266582c 6 + 749954442094c 8 

p v ' 8921808896(1 + c 2 ) 5V 

+2536970204990c 10 + 3859089592210c 12 + 5130563563118c 14 + 5335972052992c 16 + 4426939814290c 18 

+2910344048750c 20 + 1498699974850c 22 + 592140524690c 24 + 17341 7402730c 26 + 35485039590c 28 

+4530776250c 30 + 271846575c 32 ) 

271846575 , ~ Nl9 arctanc N 

+ 892T808896 (1 + C) ^^- (Al8) 
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